Adiabatic quantum pumping, magnification effects and quantum size effects of 

spin-torque in magnetic tunnel junctions 
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We study the adiabatic quantum pumping and quantum size effects of spin-torque in a magnetic 
tunnel junction within a scattering matrix approach. Quantum size effects are predicted in the 
presence of a dc bias as a function of the thickness of the normal metal layer inserted between two 
magnetic layers and of the fixed magnetic layer. In the presence of ac voltages, the results for the 
spin-torque show a peculiar magnification effect and advantages of spin-torque pumping in actual 
devices are also discussed. 

PACS numbers: 73.23.-b,72.25.Pn,75.60. Jk,72.15.Qm 



I. INTRODUCTION 

Magnetic tunnel junctions (MTJs) are layered struc- 
tures of alternating magnetic layers (generally ferro- 
magnetic) and nonmagnetic layers (insulating or nor- 
mal tunnel barriers) which have recently attracted a 
lot of attention because of magnetoresistance (MR) and 
spin-valve effects*^— MTJs based on epitaxial MgO 
barriers^ are used in the magnetic random access 
memory (MRAM) devices that work by spin-transfer 
torque. 6,7 While initially MTJs have attracted a substan- 
tial attention for their tunnel magnetoresistance (TMR) 



propertie; 



3.8-10 



more recently, the focus has shifted 



to spin-transfer torque and current-induced magnetiza- 
tion switching.— ~— From the theoretical point of view, 
spin-transfer properties have been studied extensively 
in spin valve structures based on various models- 
while for tunneling structures such studies are still rela- 
tively few i 15 ' 17 " — The spin-transfer torque has been an- 
alyzed in MTJs by first-principles electronic structure 
calculations 3 -! or by Boltzmann equation 3 ^ and recently 
it has been revisited in the Stoner model by scattering 
theory^ 3 - and in spin-valves by first principles with the 
aid of scattering wavefunction 3 ^.. Ab initio studies of the 
spin toraue in metallic GMR junctions are also reported 
in Ref. [351] . Despite the large amounts of spin-torque re- 
lated papers, among the aspects which require further 
attention are quantum size effects due to thin normal 
insertions layers and alternative ways of generating mag- 
netic torques. 

In this paper we will analyze the spin-transfer torque 
in a quasi-one-dimensional magnetic tunnel junction in 
which two magnetic regions (F1/F2) are separated by 
normal nonmagnetic spacers (NM). We will focus on a 
quantum pumping mechanism as a nonconventional way 
of generating spin-torque and the advantages of such 
mean compared to the conventional generation by ex- 
ternal dc bias will be discussed. In particular, a magni- 
fication effect of the spin-torque will be predicted in the 
presence of ac bias. In the case of a conventional genera- 
tion of spin-torque by an external dc bias, quantum size 



effects in the spin-transfer torque on a free layer will be 
analyzed and how to get information on the polarization 
at the interface will be illustrated. 

For our analysis we choose a fully quantum mechanical 
treatment of transport based on the scattering approach 
in a ballistic regime. The motivation for concentrating 
on the ballistic regime, i.e. on structures in which the 
transverse wavevector is conserved during transport, de- 
rives from the evidence of quantum oscillations observed 
in various FNIF structures. The ballistic regime is also 
characterized by a spin-diffusion length l s and a mean- 
free-path l m larger than the whole microstructure. 

The organization of the paper is the following: In Seclrfl 
we introduce the model Hamiltonian and present the 
scattering matrix approach generalized for the calcula- 
tion of the spin-torque in a MTJ. We then derive the 
expression of the torque components caused by an exter- 
nal dc voltage bias and the spin torque pumped via two 
ac voltages. In Sec lIIII we present the results of the spin- 
torque for the structure shown in FigJTJ Compared to 
Ref.fll] we focus here on quantum effects related to the 
finite width of the magnetic layer and on magnification 
effects of spin-torque by quantum pumping. 



II. 



THE MODEL AND FORMALISM 



Our system is shown in FigJTJ it is a multilayer struc- 
ture connected to two external leads in which two mag- 
netic layers are separated by nonmagnetic normal (NM) 
insertions (i.e. a NM/F1/NM/F2/NM microstructure). 
For simplicity, one of the layers Fl is taken to have a 
width less than the De Broglie wavelength and thus acts 
like delta barrier spin-dependent potential. The system 
Hamiltonian is the following: 

h 2 

H = -—dl + -/ 1 (x)n 1 -a + l2 (x)n 2 -a + V{x), (1) 
Am 

where ^i=\,i(x) — S ^-Bi(x) is the exchange splitting, 
rii is the unit vector in the direction of the exchange 
splitting and a = (o~ x , a v , cr z ) are the Pauli matrices. 
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In our specific model B\{x) — B\l8(x), while B2(x) = 
B-iX^Ai^i + ds,x) where the function \ is defined as 
X{xi,Xf,x) = 0(x — Xi)9{xf — x) and 9 is the Heaviside 
step function. The spin independent scattering poten- 
tial V(x), which can be controlled by means of the gates 
G1/G2 is given by: 

V{x) = X {-di,0,x)Vi + x(0,d2,x)V 2 + V t5(x). (2) 

The parameter I which multiplies the local d(x) poten- 
tials in Bx{x) and V(x) accounts for the layer finite size 
effects and is approximately equal to the size of the layer 
assumed much smaller than the Fermi wavelength 39 . 
We are interested here in the calculation of the spin- 
torque experienced by the layer PI at x = 0. 

The spin torque T is defined as time derivative of the 
electron spin, represented by the operator s. This yields 
the total spin torque as T = —(i/h){s,H} = Ti + %, 
where the Tj — 2(jj/h)rij x s, j = 1,2, is the torque 
experienced by the j-th FM layer. However, when the 
magnetization direction of a given layer is fixed (fixed 
layer) only the magnetization direction of the other, the 
so-called free-layer, can be affected by the local torque in- 
duced by the gradient of a spin polarized current. Thus 
in the following we focus on the spin-torque on the free 
layer. This torque can be measured by tunnel magne- 
toresistance experiment. 

In our set-up of FigfTJ the layer Fl represents the free- 
layer, being F2 the fixed layer whose magnetization di- 
rection is hi = (0,0,1). The torque generated on Fl, 
i.e. 71, lies on the plane perpendicular to the magneti- 
zation direction hi — (sin(0), 0, cos(0)) of the free-layer 
since hi ■ 71 =0. The projection of 71 parallel and per- 
pendicular to the free-layer can be expressed in terms of 
the following set of basis vectors: 



v± = 



722 x Til 



\n 2 x ni| 
hi x (fi2 x hi) 
\hi x (h 2 x hi)\ 



(3) 



-x cos(#) + z sin(0), 



where x, y, z are unit vectors along the direction of the 
cartesian axis. The torque acting on the free-layer can 

thus be decomposed as T\ — TiV\ \ + Ti~v±, where 



Ti 1 



Tl-V\\=Tl,y (4) 

Ti'0± = -Ti, x cos(6>) + Ti, z sin(0). 



To calculate the torque components acting on the free- 
layer at x = we make use of the scattering matrix 
approach and of the following definition: 



(Ti^x = 0)) = ^ h(x = 0) x (s(x = 0)) 



(5) 



where /i = {x, y, z}, h(x = 0) = hi, while (s(x = 0)) rep- 
resents the quantum average of the spin density operator 
on the free layer. 



This quantum average can be evaluated by using the 
expression of the electron field operator within the scat- 
tering approach of Ref . [4(| : 



V a (x,t) = ^2 / dE Pa (E)exp 



E 



\a)x (6) 



Akx n a 



-ikx m 



where p a (E) = [y / 2nhv a (E)] 1 is the density of states 
of the external lead a = 1,2, while v a (E) is the velocity 
of the electrons with wave vector k{E) — \/2mE/h. The 
scattering operators a Q , b a for the incoming and outgoing 
states, respectively, are related by the scattering matrix 
S through the relation b a = ^« S^a 13 and the notation 
a a ,b a stands for the following spinorial representation: 



(7) 



while a"t — 



(a? , ) , and similarly for b. Using 



the relation (a^ (E)a^,(E')) = S af3 5 a ^5(E - E')fp(E), 
f 13(E) being the Fermi function of the lead j3, the 
spin density can be computed as the quantum average 
(^a(x, t)%B^ a (x, t)). Since we are interested in the spin 
density on the free layer only the contribution from the 
closest lead, the left one (a = 1), can be retained. Using 
this expression the spin torque can be evaluated via 
Eqs.fll-©. 

Two different ways of generation of spin torque can be 
considered: one relies on the conventional application of 
dc external voltage bias V, the other is the adiabatic 
quantum pumping mechanism. The corresponding ana- 
lytic expressions will be derived below. 
Let us finally comment on a physical difference between 
the in-plane and out-of-plane spin torque components. 
As shown they originate from different spin vector com- 
ponents of the spin current implying a further qualitative 
differences between them. A natural difference is that 
the out-of-plane component is present in equilibrium, i.e. 
at zero bias, when the external leads present a spin de- 
pendent energy spectrum. This difference can also be 
understood on general symmetry grounds ( see Ref. [ill] ) . 



A. Spin torque by dc voltages 

The /i-th component of the spin density (s^) in the 
approximation of a constant density of states at the Fermi 
level, p(E) r-j p(Ep), can be written as^: 



1 



Attvf 



E 



dETr{S 1 ^(E)a^(E)}f p (E), (8) 



vf being the Fermi velocity. When a dc external volt- 
age V is applied to the microstructure fp (E) —> fp(E± 
eV/2), and Eq.© can now be used to compute the torque 
components acting on the free-layer. To linear order in 



3 



V one gets: 

7i x = 



Z7T 



(9) 



wT 



Tr[(sm(6)<j z - cos(6)a x )S 12 S* 12t ] 



where w = eV/2, while the dimensionless parameter T 
is r = Recalling that 7 = gfj, B B x £/2, T can 

be rewritten as the Zeeman energy (normalized to the 
Fermi energy Ep) of the free- layer rescaled by the nor- 
malized effective length k F £ (i.e. V = (k F £)^§^). 
Apart the spin torque, one can define the torkance as 
the linear response to a small variation of the external 
bias Sw = eSV/2, (dT l _Jdw)Sw, and its expression is then 
T P =VW2). 



B. Quantum pumping of spin torque by ac 
external gates 
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FIG. 1: Representation of the NM/F1/NM/F2/NM system 
and of the respective potential energy. Fhe spin current flows 
along the a;-direction, while the magnetizations hi and fi2 
belong to the x — z plane. 



Quantum pumping 3 - is a well known quantum effect 
for charges. In a charge quantum pump a dc parti- 
cle current is generated by the ac adiabatic modulation 
of at least two out-of-phase independent parameters of 
the system in absence of bias. In our calculation we 
will use the idea that in a magnetic layered structure 
a pumping mechanism can generate spin currents other 
than charge currents, and thus a spin-torque is generated 
on a magnetic layer by the gradient of spin current, or 
equivalently by the local spin-density (see Eq.flSJ). Re- 
cently quantum pumping has been proposed as an addi- 
tional control of the spin flux in absence of external dc 
bias^l. Focusing on the microstructure of FigfT] and ap- 
plying the idea of pumping, we modulate harmonically 
in time the barriers heights by two top gates GI and 
G2. When the gates are varied adiabatically in time the 
scattering matrix depends on time via the two varying 
external parameters as S{t) = S{X\{t), X 2 (t)), where 
Xi(t) = X° + X? sm{uit + <fi) (i = 1,2). In the weak 
pumping regime, i.e. when X£ 2 ^ -^121 the scattering 
matrix can be expanded as follows: 



S a0 (t) ~ Sf + Y Sfe ir > u 
i;=±i 



(10) 



where to = 2nv is the pumping frequency and the matri- 
ces S'™ /3 are given by 



Sn 



X?(d Xl S) +X%e ir >«>(d X2 S) 



(11) 



if being ip 2 — ipi. The Fourier transform of (110[) is then: 
S^(E) = 2n[SfS(E) + ]T SfS(E + rp)]. (12) 

»7=±1 

Using this equation in evaluating the spin density, its 
components per unit of area are given by: 

4 i^E / dETriSW^s^ME+nu). (13) 



Since no bias is present between the leads, i.e. fi(E) = 
/2(E) — f(E), and in the zero temperature limit, the 
/i-th component of spin-density per unit area to leading 
order in the adiabatic frequency uj is: 



4:77V F 



(14) 



H>1 



Substituting (fT4f in ([5]) we obtain the explicit expressions 
of the pumped torque components acting on the free- 
layer: 



V 1 



hioT 
~8t7 

8tT 



X?X% sin(^) Y Tr[Af + if] (15) 



X?X$ sin(^) Y, Tr[sm{6){Af + A™) 



s{9){Af + Al^)] 



where the quantity A<f = i(d X2 S al3 )la tl {dx 1 S al3 )o has 
been introduced. 

The equivalent of the torkance in the dc case is obtained 
for the pumping case by T M = T^/{hv/2). In the fol- 
lowing we introduce the dimensionless potential barriers 
n = (i = 1,2) and r = (kp£)Vo/Ep, the normal- 
ized Zeeman energy of the fixed layer h — 9t 2 L E B2 and the 
dimensionless distances kpdi. 



III. RESULTS 

In the following we present the results for the torkance 
components per unit of area, or equivalently of the spin 
torque normalized by eU/2 in the dc case and by hv/2 in 
the pumping case. 
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A. dc case 

In FigfJ] the torkance components (TjijJ are plotted 
in units of area as a function of the Zeeman energy h of 
the fixed layer for the remaining parameters: T = 0.5, 
r% = ri = ro = 0, 9 = tt/2, kpd\ = kpdi = 3 and 
kpds = 1. In absence of scattering potentials along 
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FIG. 2: Torkance T± (full line), 7j (dashed line) plotted as 
a function of the Zeeman energy of the fixed layer h. The 
remaining parameters have been fixed as follows: T = 0.5, 
n = r'2 — ro = 0, 6 — 7r/2, fc^rfi = fcFafe = 3 and kpdz = 1. 

the transport direction (i.e. r% = r2 = J"o = 0) and in 
presence of the exchange interaction only, the parallel 
component of torque satisfies the relation Xj| > Tj_, as in 
conventional spin-valves. The component Tj_ becomes 
relevant for r^ ^ (not shown here). Furthermore, 
if the fixed layer is made of a weak ferromagnet (i.e. 
h < 0.2) the torkance components present a linear 
dependence with respect to h, while deviations from the 
linear behavior are observed for increasing values of the 
Zeeman energy. In particular, for a critical value of the 
Zeeman interaction, close to h ~ 0.75, the perpendicular 
component of the torkance T± is totally suppressed, 
while Tj| is the only relevant component. 

In Fig(3]we plot the torkance components (T||,j_) as a 
function of the width of the fixed layer fcpGfa for the re- 
maining parameters: T — 0.5, r\ = r-i = ro = 0, = tt/2, 
kF di = kpd2 = 3 and h = 0.2. In agreement with what 
found above, in absence of scattering potentials along the 
x direction the parallel component of the spin torque Tj | is 
larger than T± over a large range of the fixed layer width. 
As shown in the figure, the torkance presents a charac- 
teristic oscillatory behavior as a function of the width of 
the magnetic layer. These oscillations can be regarded as 
a quantum-size effect. They reflect the perfect ballistic 
regime of electron transport through the multilayer struc- 
ture. The physical mechanism behind the oscillations is 



the interference effect of the electrons propagating across 
the non- magnetic/magnetic interface from the left lead 
to the right lead and electrons propagating backwards 43 . 
The particular value of the oscillation follows from the 
values of the spin-dependent Fermi wavevector. Fur- 
thermore, the behavior found in Figl3]is similar to that 
found in Ref.— (see Figs. 2, 3 of the cited work) where the 
torkance of a spin-valve was analyzed by ab initio calcu- 
lation. In our case, differently from [3l|, we observe large 
oscillations of the Tj | component around the mean value 
due to the quasi-one-dimensional character of our sys- 
tem. In the two-dimensional case, not considered here, 
we expect sizable changes. In fact, when the integration 
over the Fermi surface is performed taking into account 
all the incident directions of the electrons momentum a 
reduction of the amplitudes of the oscillating part of the 
torkance is expected. Very remarkably, the oscillatory 
behavior shown in Figj3] and particularly evident for the 
T| | component displays slow and fast scales of oscillation 
with respect to the fixed layer width. This can be seen 
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FIG. 3: Torkance T± (full line), Tj| (dashed line) plotted as 
a function of the normalized width kpd3 of the fixed layer 
F2. The remaining parameters have been fixed as: T = 0.5, 
ri = r-2 = ro = 0, 9 — tt/2, kpdi = kpd^ = 3 and h = 0.2. 

by looking at the behavior of Tn as a function of the fixed 
layer width kpds for different values of the Zeeman en- 
ergy h and fixing the remaining parameters as in Fig|3] 
The results are shown in FigfJ] where Tn has been plot- 
ted vs kpds for values of h ranging from h = 0.15 (lower 
curve) up to h = 0.45 (top curve) with a step of 0.05. The 
analysis of the lower curve in Fig|4] shows an oscillating 
behavior vs kpds characterized by fast oscillations with 
a frequency f2+ , the amplitude of the signal being mod- 
ulated by a curve with frequency f2_. As shown in FigfJ] 
by increasing the Zeeman energy h the larger frequency 
f2_i_ remains almost unchanged while the smaller one 
increases. This trend is represented by the dashed line 
in FigU] . The oscillating behavior of T\\ as a function of 
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FIG. 4: Torkance Tj| plotted as a function of the normal- 
ized width k F d 3 of the fixed layer F2. The lower curve is 
plotted for h = 0.15, while the artificial h-dependent upward 
shift T u (k F d 3 ) T u (k F d 3 ) + 0.02n, where n = 0,1,... for 
h — 0.15, 0.2, ... has been given to the curves to facilitate the 
comparison. The remaining parameters have been fixed as in 
FigE 



k F d 3 can be fitted by a non-linear regression with trial 
function: 

T\\{k F d a ) = T + cos(il-k F d 3 )[Asm(Q + k F d 3 ) +(16) 
+ B cos(il + k F d 3 )}, 

To, A, B, n± being fitting parameters. In Figf5]we report 
the fitting analysis of fi_(/i) (full circles) and for com- 
parison we also plot the line CI -(h) = h. The analysis of 
Figf5] shows that the slow frequency 51_ is controlled by 
the Zeeman energy h of the fixed layer and thus indirectly 
can give a measure of the polarization of the electrons 
belonging to the fixed layer. Indeed, remembering that 
the matching conditions on the electron wavefunctions 
needed to compute the S-matrix involve oscillating func- 
tions of the form exp(±ikcrd 3 ), where k a = kpy/1 + ah, 
one expects that the scattering matrix elements to the 
lowest order can be approximated via a linear com- 
bination of terms sin(fc CT d3) and cos(k a d 3 ), or equiva- 
lently by harmonic functions of argument (fa\- ± fcj,)/2. 
Since the spin torque depends roughly on the scatter- 
ing matrix elements squared, one expects that the os- 
cillation frequencies of Tim vs k F d 3 are of the form 
fl± = vT+H ± y/l — h. Assuming this relation for 
f2+, ft- one reproduces exactly the result in Fig|5]in the 
limit of small h. Moreover by considering the limit h — » 
one expects fi_ — > while only the oscillations with fre- 
quency f2+ survive. Indeed this is found in Fig|6] where 
the torque components as a function of non-magnetic 
layer width k F d2 is shown. Indeed, the figure presents 
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FIG. 5: fi- as a function of h (full circles) obtained by 
non-linear fit procedure of the curves shown in Fig(4] The 
dashed line, inserted for comparison, represents the linear 
model Q-(h) = h. The full line interpolates the computed 
points and shows the deviation from a simple linear model. 



an oscillating pattern of period 2tt/CI + = ir, the modula- 
tion with smaller frequency being totally absent. 
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FIG. 6: Torkance T± (full line), Tj| (dashed line) plotted as a 
function of the normalized spacer width k F d2- The remaining 
parameters have been fixed as follows: T — 0.5, n = r2 = 
r = 0, 6 = tt/2, = 4, fc F rf3 = 2, h = 0.25, y> = vr/2. 

In FigJ7]we show the torkance components as a func- 
tion of the angle 9 between the magnetizations of the 
regions Fl and F2. Apart from the sinusoidal behavior 
of the curves with respect to 9, the presence of scatter- 
ing potentials along the transport direction (lower panel) 
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FIG. 7: Torkance T± (full line), Tj| (dashed line) plotted as 
a function of the angle 6. The remaining parameters have 
been fixed as follows: F — 0.5, k F di = k F d 2 = 3, = 10, 

h = 0.45, tp — 7r/2 and ri = r2 = ro = (upper panel), or 
ri = 0.1, T2 = 0.05, ro = 0.2 (lower panel). 



enhances the perpendicular component of the torkance if 
compared with the one obtained in the absence of scat- 
tering potentials (upper panel). Furthermore the vs 
curves present an in-phase behavior which can be al- 
tered by changing the relative size of the microstructure 
due to the oscillating nature of the torkance as a function 

of kpd%. 



B. Pumping case 

Here we analyze the spin torque generated on free layer 
by a quantum pumping mechanism. In particular, we 
adiabatically modulate in time the electrostatic poten- 
tials ri{t) as: 



where i — 1,2, and make the following choice for the 
pumping phase <px = 0, ifi = tt/2. 

In Fig. [3] we plot the torkance T]|_j_ as a function of the 
fixed layer width kpd% fixing the remaining parameters 
as: T = 0.5, r\ = r% = r = 0, = r% = 0.1, 6 = tt/2, 
kpdi = kpdi = 3, h = 0.2, ip = tt/2. By comparing FigJS] 
with the analogous figure obtained for the dc case (i.e. 
FigOJ) one observes that the pumping procedure modi- 
fies the sign of the mean values of torque components 
and simultaneously enhances the relative amplitudes of 
the oscillating patterns. As discussed above, in presence 
of scattering potentials along the transport direction, we 
expect the perpendicular component of the spin torque 
to be more relevant and, compared to the dc case, this 
is particularly true for the pumping case where two scat- 
tering potentials are modulated in time. Again one can 
distinguish two frequencies of oscillation of the torque 
components Q± whose values are the same obtained in 
the dc case. Apart from the changes described above, 
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FIG. 8: Torkance T ± (full line), T\\ (dashed line) plotted as 
a function of the normalized fixed layer width kpd?,. The 
remaining parameters have been fixed as follows: T — 0.5, 
r? = r% = r = 0, r^ = r£ = 0.1, 6 = tt/2, k F d 1 = k F d 2 = 3, 
h = 0.2, (p = n/2. 

the quantum pumping mechanism may induce relevant 
effects on the spin torque experienced by the free layer 
due to the parametric derivatives of the scattering ma- 
trix that appear in (p~5|) . This is clearly seen in Figj9] 
where the torque components Tji j_ are plotted as a func- 
tion of the metallic spacer width kpdi- Differently from 
the dc case shown in Fig[51 in FigJS] we observe a lin- 
ear increasing of the oscillation amplitudes of the torque 
components and the lost of periodicity of vs kpdi. 
This is a magnification effect of the torque pumped in the 
system that can be exploited in current experiments. Let 
us note that the value kpdi = 20 corresponds to a spacer 
width of ~ 3.18Ai? thus our balistic treatment is still ap- 
propriate. The linear increasing of the oscillation ampli- 
tude can be naively explained observing that the pumped 
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FIG. 9: Torkance T x (full line), Tj| (dashed line) plotted as a 
function of the normalized spacer width fc F d2. The remaining 
parameters have been fixed as follows: Y = 0.5, r\ — r% = 
r = 0, jf = r% 
h = 0.25, ip = tt/2 
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FIG. 10: Torkance T± (full line), T N (dashed line) plotted as 
a function of ro. The remaining parameters have been fixed 
as follows: F = 0.5, r? = r£ = 0, rf = r£ = 0.1, = tt/2, 
fc F di = 4, k F d 2 = 10, fc F d 3 = 2,h = 0.25, </? = tt/2. 



torque is related to the derivative of a periodic function 
of fcFC?2, F — asin(vl — r-ikpd-i) + 6cos(\/l — T2kpd2) 
with respect to r2 (barrier height) that gives a coeffi- 
cient kpd,2- Concerning the dependence of the torque 
components on the scattering potential amplitude along 
the transport direction, we can observe in Fig llOl an in- 
creasing of the perpendicular component of spin torque 
Tj_ as a function of the barrier located on the free-layer 



ro, while the component Xji becomes very small. When 
the transparency at the free-layer ro becomes small (i.e. 
for high values of ro) both the components of the spin 
torque decrease due to a suppression of the spin-fluxes. 
In Fig llll we present the torque components as a function 
of the Zeeman energy h of the fixed layer for fcj?c?3 = 6 
(upper panel) or kpd^ — 5 (lower panel) and fixing the 
remaining parameters as follows: T = 0.5, r® — r% = 0, 
r = 2, rf = r'i = 0.1, 6 = tt/2, k F d x = 3, k F d 2 = 3, 
if = tt/2. One can notice a strong dependence on the 
width of the fixed layer and, in particular, one observes 
a sign reversal of T±_ at varying h (lower panel). This 
behavior can be understood by the analytical expression 
of the torque which is an oscillating function of argu- 
ment hkpd^. As shown above a certain value of h, i.e. 
h c ~ 0.75, the magnetic barrier height increases reduc- 
ing the transmission of spins across the layer and thus a 
suppression of the spin-torque is observed (over-barrier 
reflection). 
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h 

FIG. 11: Torkance T± (full line), (dashed line) plotted as a 
function of the normalized Zeeman energy h. The remaining 
parameters have been fixed as follows: F — 0.5, r? = r§ = 0, 
ro — 2, r" = r2 = 0.1, 9 — tt/2, fc F di = 3, fc F c?2 = 3, 
ip — tt/2, while the width of the fixed layer is set to fc F d3 = 6 
(upper panel) or fc F d3 = 5 (lower panel). 
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FIG. 12: Torkance T± (full line), Tj| (dashed line) plotted as 
a function of the angle 0. The remaining parameters have 

kpdi = fcFt?2 = 3, kpds — 10, 



r° 2 







been fixed as follows: F — 0.5 
h = 0.45, p = tt/2, rf = r£ = 0.3 and r? 
(upper panel); F = 0.65, fc^cii = 3, kpd.2 = 15, fcj?d3 = 10, 
h = 0.85, if = tt/2, = = 0.3 and r? = r% = r = 
(lower panel). 



In Fig[T2]we report the torkance as a function of the 
angle between the two magnetizations h\ and fi2, while 
setting the other parameters as in the figure label. As 
shown, the perpendicular component of the torque close 
to 9 = tt/2 is ~ 30% of the parallel component due to 
the scattering potentials r^ 2 , while a sinusoidal behav- 
ior with respect to 6 is observed for both curves. Again 
one observes a magnification effect of the spin torque by 
increasing the width kpdi of the non- magnetic metallic 
layer between the two ferromagnetic layers Fl and F2 
(the values of the torque increases by a factor 10 2 ). The 
magnification effect reported above can be easily under- 
stood as a combined effect originated by the pumping- 
induced magnification related to the increasing of the size 
fcp c?2 of the interstitial layer further enhanced by exploit- 
ing strong polarized magnetic layers F1/F2 (i.e. by set- 
ting h — 0.85 and T = 0.65 instead of the values h = 0.45, 
r = 0.5 used in the upper panel of Fig fTS")) . It is worth 
to mention that the torque per unit of area induced by 
the quantum pumping within the weak pumping regime 



analyzed here is typically a fraction in unit of /leV/area. 
At pumping frequency of vq = 300 MHz, the energy scale 
normalizing the torkance Hvq/2 is given by 0.1 [J.eV, i.e 
the ~ 12.5% of the maximum value (~ 0.8 fieV / area) 
obtained within the dc case in Ref.[44j by considering 
a MTJ. Thus in our simulations, considering vq = 300 
MHz, we have verified that the maximum value of the 
torque obtained by a weak pumping procedure ranges 
from 3% up to 6% of the maximum value of the spin 
torque obtained in Ref.j44|. Higher values of the spin 
torque produced by a quantum pumping mechanism can 
be obtained beyond the weak pumping regime. 



IV. CONCLUSIONS 

We analyzed the spin-torque in a NM/F1/NM/F2/NM 
microstructure by a scattering matrix approach and con- 
sidering different mechanisms for its generation: 1) the 
spin torque induced by a dc voltage V applied to the 
whole system; 2) the spin torque activated by the quan- 
tum pumping technique. While the first method is widely 
studied in the present literature, the second one based 
on the quantum pumping has been initially proposed in 
Ref.[38j and a complete analysis has been performed in 
the present work. In particular, we have analyzed the 
quantum size effects induced by the finite width of the 
ferromagnetic fixed layer F2 both when the system is 
forced by a dc bias and in the pumping case and found in- 
teresting features related to the quantum pumping mech- 
anism. In the case of an external dc bias, the most evi- 
dent feature is the presence of an oscillatory behavior (de- 
tectable also in the pumping case) of the Tii ± vs kpd^ 
curves characterized by frequencies fl± directly related 
to the Zeeman energy of F2. These oscillations reflect 
the perfect ballistic regime of electron transport across 
the whole system and can give important information on 
the polarization at the interface of the magnetic layers. 
Indeed, the spin torque arises either as an interference 
effect between spin up electrons propagating across the 
ferromagnetic region from the R lead to the L lead and 
spin-down electrons propagating backwards or have to be 
ascribed to quantum well states, i.e. to an interference 
effect in a single spin channel. The finite layer width ef- 
fects described within the dc case by using our theory 
present qualitative agreement with recent studies on the 
spin torque generated in a Cu/Fe/MgO/Fe/Cu tunnel 
junctional and in a Cu/Co/Cu/Ni/Cu system^. Then, 
we have proposed a parametric quantum pumping of spin 
torque. The pump works by means of two external gates 
able to produce out-of-phase voltage modulations on two 
nonmagnetic regions attached to the free layer (i.e. a 
thin ferromagnetic region Fl). The underlying idea is 
that in a magnetic layered structure a pumping mech- 
anism can activate spin currents other than charge cur- 
rents and thus a spin-torque is generated on the magnetic 
layer subject to the spin-current gradient. This peculiar 
way of generating spin torque is strongly affected by the 
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dependence of scattering matrix of the microstructure 
on the pumping parameters (i.e. the external voltages 
controlled by the gates G1/G2). As a consequence of 
this parametric dependence of the scattering matrix, a 
peculiar magnification effect of the perpendicular com- 
ponent of torque has been predicted, the latter feature 
being particularly appealing to test the proposed theory 
using MTJs or exploitin g a modified system similar to 
the one described in Ref.[45|. Indeed, by increasing the 



width of the nonmagnetic spacer kpd,2, a magnification 
of the torque components has been observed, differently 
from the dc case where the Tj| t ± vs /cf^2 curves present 
a simple oscillating behavior. Our estimate of the spin 



torque induced by the weak pumping (using a pumping 
frequency v = 300MHz) is 6 — 7% of the one obtained 
conventionally using dc voltages, nevertheless the effects 
of magnification can be efficiently exploited beyond the 
weak pumping regime to obtain values of 2]i j_ similar to 
the one observed in the dc case. Apart from the tech- 
nological motivations supporting our work, the quantum 
pumping of spin torque can be considered as the proto- 
type of a new class of quantum pumps able to pump a 
vector (i.e. the torque) instead of a scalar (i.e. the elec- 
tron/hole charge) and can be relevant to further test the 
quantum effects in nano-electronics. 
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